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ABSTRACT 

We estimate that the moment of inertia of star A in the recently discovered double pulsar 
system PSR J0737-3039 may be determined after a few years of observation to something like 
10% accuracy. This would enable accurate estimates of the radius of the star and the presure of 
matter in the vicinity of 1 to 2 times the nuclear saturation density, which would in turn pro- 
vide strong constraints on the equation of state of neutron stars and the physics of their interiors. 

Subject headings: Equation of state — moment of inertia — binary pulsars stars: neutron - 
stars: pulsars 

1. INTRODUCTION 

The discovery of the double-pulsar system PSR J0737-3039 (Burgay et al. 2003; Lyne et al. 2004) provides 
physicists with a remarkable laboratory for relativistic astrophysics. Besides its implications for the rate of 
gravitational wave bursts from neutron-star coalescence (Burgay et al. 2003) and for the understanding of 
pulsar magnetospheres (Lyne et al. 2004), it could provide a measurements of spin-orbit coupling. Lyne et 
al. (2004) noted that accurate timing over a period of years could lead to a determination of the moment of 
inertia of star A. Spin-orbit coupling could be revealed either through an extra advancement of the periastron 
of the orbit above the standard post-Newtonian advance, or in the precession of the orbital plane about the 
direction of the total angular momentum of the system. Given that the masses of both stars are already 
accurately determined by observations, a measurement of the moment of inertia of even one neutron star 
could have enormous importance for neutron star physics. 

Despite the fact that over a thousand neutron stars have been discovered in radio and X-ray observations, 
and accurate masses have been determined for a dozen or so neutron stars in radio pulsar binaries (Stairs 
2004), there is relatively little observational information about their radii or their interior physics. Currently, 
one source of candidates for revealing radii are those neutron stars for which thermal x-ray and optical 
radiation have been observed (see Page et al. (2004) for a review). Nevertheless, these stars have not yet 
provided a clear radius determination, because of ambiguities due to atmospheric re-processing, interstellar 
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absorption, and distances. Moreover, the inferred radius actually refers to the radiation radius, = 
R/y/l — 2GM/Rc 2 , which is mass-dependent; unless independent mass measurements of the same stars 
become available, such as would be provided by a redshift, the radii themselves remain unknown. However, 
most of these sources have featureless spectra and no secure redshifts have been obtained to date from them. 

Another possibility are X-ray bursts from the surfaces of neutron stars. Recently, two lines observed 
in an X-ray burst spectrum have been suggested to be H- and He-like Fe lines and imply a redshift of 0.35 
(Cottam, Paerels & Mendez 2002). This inference has been given additional credibility by the detection 
(Villarrcal & Strohmayer 2004) of a 45 Hz rotational frequency for the neutron star, EXO 0748-676. A 
low spin rate is consistent with the observed equivalent widths of these lines if their identifications with Fe 
are correct. Villarreal & Strohmayer (2004) show that this consistency holds if the neutron star radius is 
in the range 9.5-15 km (corresponding to masses in the range 1.5-2.3 M Q ). Since this star is a member 
of an eclipsing binary, an independent mass measurement might be possible, as well. These techniques 
hold promise of further constraining the radii of this star, and could be extended to other X-ray bursters if 
redshifts from them can be observed. 

At present, however, the constraints on the neutron star radius from X-ray bursters and thermal radia- 
tion from cooling neutron stars are relatively weak. We demonstrate that securing a value for the moment of 
inertia for a star in a radio binary pulsar system will provide important constraints on the radius of the star 
and the equation of state of neutron-star matter. Dimensionally, the moment of inertia is proportional to 
the star's mass times its radius squared, so that a measurement of the moment of inertia to a given accuracy 
provides approximately twice the accuracy for a radius identification. 

We begin in §2 by exploring the remark in Lyne et al. (2004) that the moment of inertia can be deter- 
mined from measurements of spin-orbit coupling. Two observable effects pf spin-orbit coupling, precession 
of the orbital plane and an extra contribution to the advance of the periastron, are investigated. In §3, the 
case of PSR 0737-3039 is investigated. We conclude that a 10% measurement of the moment of inertia of 
star A in this system might be possible with observations extending over a period of a few years. In §4 
section, we show how such a measurement leads to an estimate of the radius of the star and the density of 
neutron star matter in the vicinity of the nuclear saturation density, This could have crucial implications for 
delimiting the equation of state. The conclusion section, §5, compares the PSR0737-3039 system with other 
rclativistic binaries, and possibilities for the future of this technique are discussed. 

2. OBSERVABLE SPIN-ORBIT EFFECTS 

There are two kinds of spin-coupled precession effects in a binary system of compact stars: spin-orbit 
and spin-spin couplings. Spin-orbit coupling leads to a precession of the angular momentum vector L of 
the orbital plane around the direction of the total angular momentum J of the system. This is sometimes 
called geodetic precession, and is related to the Thomas precession of atomic physics. Since the total angular 
momentum J = L + Sa + Sb is conserved (at this order), there are compensating precessions of the spins 
Sa and Sb of the two stars. Since the orbital angular momentum dominates the spin angular momenta in 
binaries, the geodetic precession amplitude is very small while the associated spin precession amplitudes are 
substantial. In addition to geodetic precession, spin-orbit coupling also manifests itself in apsidal motion 
(advance of the periastron). Spin-spin coupling is negligible in close binary systems because \L\ » \Sa\, \Sb\- 



- 3 - 



According to Barker & O'Connell (1975), the spin and orbital angular momenta evolve according to 

3. G(4M I + 3M_ t ) z g 

_ 2M 1 aV(l- e 2)3/2 ixi «- W 

~ 4" 2M ^ 3 c 2 (l-e 2 ) 3/2( |i| 2 

where the superscript SO refers to the spin-coupling contribution only (there are also first- and second-order 
post-Newtonian terms, 1PN and 2PN, respectively, unrelated to the spins, that contribute to this order). 
Here a is the semimajor axis of the effective one-body orbital problem (sum of the semi-major axes of the 
two stellar orbits), e is the eccentricity, and Mj and M_j refer to the masses of the two binary components 
(we use the notation that i = A,B and — i = B,A). To this order, one may employ the Newtonian relation 
for the orbital angular momentum: 



P - - _ ^ t^Sl (A) 



|i|= P M ~MaM b] I (3) 

where P is the orbital period and M + Ma + M B . Then, from Eq. 1, the spin precession periods are 

2c 2 aPM(l - e 2 ) 
GM_,(4JW, + 3M_ 

which arc not identical for the two components unless they are of equal mass. Note that the spin precession 
periods are independent of the spins. Also note that if the spins are parallel to L, there is, first, no spin 
precession, and second, the spin-orbit contribution to the advance of the periastron is in the sense opposite 
to the direction of motion. 

The spin precession leads to two observable effects. First, as the spin axis change orientation in space, 
the pulsar beams will sweep through changing directions in space. In many cases, this will lead to the 
periodic appearance and disappearance of the pulsar beam from the Earth. Second, since total angular 
momentum is conserved (to this order), the orbital plane will change orientation. This will be observed as 
a change in the inclination angle i. 

Damour & Schaefer (1988) have considered the question of how these effects affect the timing of binary 
pulsars. For the change in inclination, they find 

di G 7r ^/ i (4M i + 3M_ i ) . . , R . 

Jt = ^ (T-^F 1 ? M^p t sm 6i cos * ' (5) 

where 9i is the angle between Si and L and 4>i is the angle between the line of sight to pulsar i and the 
projection of Si on the orbital plane. These angles follow the convention of Jenet & Ransom (2004), but 
other references employ fa — 90 in place of fa, cf., Wex & Kopeikin (1999). Also, we used |Sj| = 2irIi/Pi 
where Ii and Pi are the moment of inertia and the spin period, respectively, of component i. If we specialize 
to the case in which the spin of one component is much greater than that of the other, \Sa\ >> \Sb\, the 
amplitude of the precession of the inclination of the orbital plane is given by the change in L needed to 
compensate changes in Sa- We define the angle between L and J to be 9 — A and the angle between Sa 
and J to be A. Since \Sa\ « \L\, one has \9 — A| << |A|. Using the fact that J ~ L + Sa, one finds 
\L\ sin(8 — A) ~ \Sa\ sin A ~ \Sa\ sm9. Thus, the amplitude of the change in the orbital inclination angle i 
due to A's precession will be 

c \S"a\ . a IaM P 

~\lT Sm " a%Jlf fl (l-^ K Sm (6) 
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This will cause a periodic departure from the expected time-of-arrival of pulses from pulsar A of amplitude 

° 1 a = ^rr _ ^ cosz = ~tj — 5-5-smvAcosi, (7) 

M c c NlAa ±A 

if one can assume the orbital eccentricity is small. 

For the advance of the periastron, the ratio of the spin-orbit and 1PN contributions is (Damour & 
Schaefer 1988) 

A^ = - 6(l- e »)V»Ma» £ MiPi (2 C ° S 9i + COt * ^ 9i ^ (8) 

In the case that \Sa\ >> \Sb\, only the i = A term contributes substantially. For comparison, the ratio of 
the 2PN to 1PN contributions is (Damour & Schaefer 1988) 

i 2PN GM ^ (\ M t fMi\ 2 ] 2 l 46M 4 10 
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where both i = A and i = B terms contribute. The second line of Eq. 9 is valid in the case that Ma = Mb- 

In §4 we demonstrate that a very useful constraint on the equation of state (EOS) can be made if the 
moment of inertia can be determined to about 10%. In practice, it is expected that the binary components 
will have approximately equal masses but differing spin periods. Therefore, the spin-orbit effects will be 
dominated by the more radidly rotating star, A. Besides being functions of the parameters Ma, Mb, P, a, e 
and i, the observables St a and A p are also functions of 9 a and (pA- Therefore, extraction of I a from 
these observables requires that additional information about the orientation of the spin of A be available. 
Fortunately, if A is observed as a pulsar, observations of the beam geometry and its precession can provide 
this information. 



3. APPLICATION TO PSR 0737-3039 

The observational parameters for the system PSR 0737-3039 are summarized in Lyne et al. (2004): 
M A 1.34 Mq, M b ^ 1.25 M Q , a/c = 2.93 s, e ~ 0.088, P A ^ 22.7 ms, P B ~ 2.77 s, and P ~ 0.102 day. 
We therefore observe that P p a — 74.9 yrs and P p b — 70.6 yrs. With these parameters, we can form the 
useful combinations 

^-=4.32HT 6 , - (7.74 lO" 11 ) I A ,so , = 3.88 10 5 (10) 

ac z Ma z Pa 

where I Am = Ia/(S0 M q km 2 ) is a typical value for the moment of inertia (see §4). 

Since Pb/Pa = 122, we can ignore contributions to the precession from pulsar B. The orientation of the 
spin axis of A relative to the orbital plane has been estimated through modeling of the intensity variations 
of pulsar B due to illumination by emission from A (Jenet & Ransom 2004). There are two solutions to 
this model: solution 1 with 6 A = 13° ± 10° (or 167° ± 10°) and 4> A = 246° ± 5°, and solution 2 with 
9a = ±90° ± 10° and 4>a = 239° ± 2°. However, solution 1 is preferred on two grounds: it is improbable 
that a supernova kick would result in a spin axis that is so strongly misaligned with the orbital angular 
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momentum, and solution 2 is also inconsistent with the misalignment angle between Sa and the magnetic 
dipolc axis, measured to be a few degrees (Demorest et al. 2004). We assume solution 1 for the following 
discussion. The orientation of B has not been estimated so far, but this is irrelevant since the spin of pulsar 
B can be ignored. 

Recent determinations of the inclination angle are i — 88.4° (Ransom et al. 2004) and i = 

90.26° ± 0.13° (Coles et al. 2004). We shall use the latter as being the more precise. The facts that the 
orbit is seen nearly edge-on and that Sa is only slightly misaligned from L makes this a special case in 
which the amplitude of the timing change produced by the orbial plane precession will be extremely small, 
StA — (0.17 ± 0.16)7,4,80 M s - Not only is the magnitude very small, but the large relative uncertainties in 
both cos? (| cosi| ~ 0.0045 ± 0.0023) and in sm0A (s'mdA — 0.22±0.17) combine to give a huge uncertainty 
in StA- Current technology allows timing accuracies at the level of tenths of microseconds, limited by the 
quality of clocks and the characterization of radio antenna signal paths (Lorimer 2001). Given this sensitivity, 
observation of the changing inclination due to orbital precession seems marginal. Even if timing sensitivities 
to order O.Ol^s were possible, errors in both i and 9a would have to be substantially reduced in order to 
measure 1a to order 10%. 

The periastron advance due to spin-orbit coupling, being proportional to sini and cos^, is, on the 
other hand, less sensitive to errors in these angles. In addition, since i ~ 90°, Eq. 8 shows that uncertainties 
in the angle 4>a become largely irrelevant. The 1PN perihelion advance is 

6tt GM 

A\pn = -, o — n~ radians per orbit, (11) 

1 — e z ac z 

or 0.294 radians per year. The periastron advance ratio is A v aI^\pn — 6-6ta6 10 _5 -^A,80- 

In practice, to measure the spin-orbit effect to 10% means subracting out the spinless pieces of the 
perihelion precession {A\pjq and A2pn) to one part in about 10 5 . Fortunately, both asini and P are well 
known, and will be even better known in the future, and the spinless periastron advance depends only on 
the total mass, not the individual masses. The total mass M = a 3 (2-7r/P) 2 is known to the accuracy of sini. 
But sini is currently known to about 2 parts in a 10 6 , using the determination of Coles et al. (2004), and its 
accuracy will also improve with time. So measurement of the spin-orbit perihelion advance seems possible. 

Nevertheless, the determination of I a from this technique will only be as accurate as our understanding 
of 9a- An uncertainty in cos9 of better than 6% has already been achieved. Further improvements in our 
knowledge of the pulse geometry, which seems likely, suggest that a moment of inertia determination accurate 
to 10%, may be possible. The epoch of periastron changes due to spin-orbit coupling by about 0.027 seconds 
per year. For comparison, the second post-Newtonian correction contributes a periastron advance time of 
about 0.02 seconds per year and must therefore be included in the pulsar timing model. 



4. CONSTRAINING THE EOS 

The estimation of a neutron star's moment of inertia from timing observations of a radio binary pulsar 
has significant implications for constraining the equation of state (EOS). In some respects, a moment of inertia 
measurement could be more useful than a radius measurement of the same accuracy. First, the neutron star 
mass in a radio binary will obviously already be known to high precision, while a radius measurement from 
observations of thermal radiation actually refers to the radiation radius R^ = R/ yjl — 2GM/Rc 2 . To obtain 
the radius and mass separately requires an estimate of the star's redshift (1 — 2GM/Rc 2 )~ 1 / 2 — 1 which so 
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far is not yet available for any thermal source. Moreover, the accuracy of radiation radius estimates from 
thermal sources are limited by uncertainties due to source atmospheric re-processing, interstellar absorption, 
and distances. In the case of X-ray bursters, the one case in which a redshift measurement may be secure, 
EXO 0748-676, no independent mass measurement is yet available, and radius information obtained from 
fitting the line profile with the observed rotation rate is relatively weak. Mass and moment of inertia 
measurements from radio binary pulsars are not sensitive to estimates of the distance, and a redshift is not 
required. Second, the range of moments of inertia for various neutron star models (including strange quark 
matter stars) is, in relative terms, larger than the predicted range of radii (see Fig. 1). This follows from 
the dimensional relation / oc MR 2 . For example, it is 30 to 240 M Q km 2 for masses larger than 1 Mq, and 
is 53 to 109 M & km 2 for a 1.4 M Q star. For comparison, the range of radii for a 1.4 M Q model is 9 to 16 km. 




M (M ) 

Fig. 1. — The moment of inertia scaled by M 3 / 2 as a function of stellar mass M. Equation of state labels are 
described in Lattimer & Prakash (2001). The shaded band illustrates a 10% error on a hypothetical I/M 3 ^ 2 
measurement of 50 km 2 Mq 1 ' 2 ; the error bar shows the specific case in which the mass is 1.34 M Q . The 
dashed curve labelled "Crab" is the lower limit derived by Bejger & Haensel (2003) for the Crab pulsar. 

Spanning the same set of equations of state as in the compendium of Lattimer & Prakash (2001), 
moments of inertia for normal neutron and strange quark matter stars are displayed in Fig. 1. The moments 
of inertia have been scaled by a factor M 3 / 2 to reduce the range of the ordinate. For most masses, the range 
in / is approximately a factor of 2 to 3. The significance of a measurement of / with ±10% accuracy is 
illustrated by the shaded band centered on the hypothetical measurement, here taken to be 7/M 3 / 2 = (50±5) 
km 2 Mq 1 ' 2 , together with an error bar located at a precisely measured mass, taken to be 1.34 M Q . It is 
clear that relatively few equations of state would survive these constraints. Those families of models lying 
close to the measured values would have their parameters limited correspondingly. 

If the equation of state does not have a large degree of softening at supernuclear densities, possibly 
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Fig. 2. — The moment of inertia as a function of the relativity parameter M/R. The curves labelled TIV, 
TVII, NIV, Inc and Buch refer to analytic solutions of Einstein's equations, while other curves refer to a 
variety of equations of state (see Lattimer & Prakash (2001) for details). The shaded band illustrates the 
relation Eq. 12. 

introduced by hyperons, Bose condensates or self-bound strange quark matter, a moment of inertia deter- 
mination furthermore permits one to estimate the neutron star radius to a relative uncertainty smaller than 
the relative uncertainty in the moment of inertia measurement. Fig. 2 shows the moment of inertia as a 
function of the relativity parameter M/R for the same equations of state displayed in Fig. 1. Several analytic 
solutions of Einstein's equations that are applicable either to normal or self-bound stars, are displayed as 
well. These solutions are all scale- free and are functions of M/R alone; hence, they cannot be displayed in 
Fig. 1. Unless the equation of state has an appreciable degree of softening, usually indicated by a maximum 
mass of order 1.6 M Q or less, or unless it is self-bound, there appears to be a relatively unique relation 
between I /MR 2 and M/R. For M/R values greater than approximately 0.07 M Q km -1 , i.e., for M > 1.0 
Mq, this relation can be approximated by 

I ~ (0.237 ±0.008) MR 2 

An analogous fit has also been suggested by Bejger & Haensel (2002): 

9 

for M/R > 0.1 M /km; however, this fit underestimates / for the largest neutron star masses. 

Simultaneous mass and moment of inertia measurements could therefore usefully constrain the radius. 
The relevant radius relation can be determined by inversion of Eq. 12. Fig. 3 shows how the radius could 
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Fig. 3. — Radius limits imposed by simultaneous moment of inertia and mass measurements, established 
from Eq. 12. Moment of inertia error bands include measurement uncertainties of 10% and systematic 
uncertainties from Eq. 12 and are labelled in units of Mq km 2 . The horizontal error bar illustrates the 
hypothetical case in which M and I are measured to be 1.34 M© and 80 ± 8 M Q km 2 , respectively. 

be constrained for selected moment of inertia measurements having 10% uncertainty. For a 1.4 M© star, 
this typically results in a radius estimate with about 6-7% uncertainty. Even in the event of significant 
softening of the equation of state, the uncertainty of the estimated radius would be degraded by no more 
than a factor of two. Of course, the accumulation of more than a single I — M pair would significantly 
enhance the constraints. 

The importance of a radius determination is that it immediately translates into a measure of the neutron 
star matter pressure near the nuclear saturation density (Lattimer & Prakash 2001) In particular, the relation 
between these quantities is of the form of a power law 

R M P(ny 1/4 = C(n,M). (14) 

Here Rm is the neutron star radius in km for the mass M and P(n) is the matter pressure in MeV frn -3 
evaluated at the density n. The constant C is parametrized by both n and M, but its dependence upon M 
is weak (approximately scaling as M _1 ' 8 ). In the case of M = 1.4 M Q , a least squares fit to approximately 
30 equations of state yielded C(n, 1.4) = 9.30 ± 0.60, 7.00 ± 0.31 and 5.72 ± 0.25 for the cases n/n s = 1, 1.5 
and 2, respectively. This relation could be made more precise by adjusting the exponent of P (Steiner 
et al. 2004), but we choose not to do so here. The pressure of neutron star matter at these densities is 
primarily a function of the density dependence of the nuclear symmetry energy. In general, therefore, we 
expect that measurement of the moment of inertia would provide estimates of pressures. Fig. 4 illustrates 
the situation for some representative equations of state for the case of M = 1.34 M Q . It is observed that the 
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Fig. 4. — Pressures at the densities n s (hollow triangles), 1.5n s (hollow circles) and 2n s (filled circles) as 
a function of the moment of inertia for 1.34 M Q stars. Solid curves show the relations between P and / 
derived by combining Eqs. (12) and (14). Dashed curves include the errors in these fits. Equations of states 
and their labels are described in Lattimer & Prakash (2001). 



phenomenological fits Eqs. (12) and (14) adequately describe all but the softest equations of state (e.g., PS 
and PAL6) for each density. The uncertainties in the estimated pressure obtained from a measured value of 
the moment of inertia are moderate, amounting to a factor of about 2 when a 10% uncertainty in the moment 
of inertia measurement is included. Nevertheless, given the fact that present estimates of the pressure of 
matter at the nuclear saturation density span a range of a factor of 6, this information will be valuable. 



5. DISCUSSION AND CONCLUSIONS 

The precession of the spins and orbital plane of the pulsar binary PSR J0737-3039 occurs with a period 
of about 75 years. Both the inclination of the orbit along the line of sight and the position of the periastron 
change due to spin-orbit coupling. The edge-on nature of this binary, coupled with the slight misalignment of 
Sa and L, probably precludes observation of the changing inclination angle, but the spin-orbit contribution 
to the advance of the perihelion, which amounts to a timing residual of nearly 0.03 seconds per year, should 
be measureable within a few years. In this case, the edge-on orbital plane and near alignment of Sa and L 
work in favor of a precision measurement of I a ■ A measure of I a to about 5-10% accuracy seems possible, 
and significant constraints on the equation of state would be forthcoming. It would also lead to more reliable 
estimates of neutron star radii and matter pressures near the nuclear saturation density than are currently 
available. 
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Tabic 1. Comparison of Binary Pulsars 





PSR 0707-3039 


PSR 1913+16 


PSR 1534+12 


References 


a, b, c 


d, e 


f, g, h 


a/c (s) 


2.93 


6.38 


7.62 


P(h) 


2.45 


7.75 


10.1 


e 


0.088 


0.617 


0.274 


M A (M ) 


1.34 


1.441 


1.333 


M B (M ) 


1.25 


1.387 


1.345 


T GW (Myr) 


85 


245 


2250 


i 


90.26 ±0.13° 


47.2° 


77.2° 


Pa (ms) 


22.7 


59 


37.9 




13° ± 10° 


21.1° 


25.0° ±3.8° 




246° ± 5° 


9.7° 


290° ± 20° 


P p a (yr) 


74.9 


297.2 


700 


5t a /lA,80 (m s ) 


0.17 ±0.16 


7.5 


5.2 ±0.8 


(^4lPJV-fA,80) 


6.6to;e x ltr5 


1.0 x 10" 5 


1.1 x 10" 5 


^Pn/^-IPN 


5.1 x 10~ 5 


4.7 x 10" 5 


2.3 x 10" 5 



a: Lyne et al. (2004); b: Solution 1, Jenet & Ransom (2004); c: Coles et al. (2004); 
d: Weisberg & Taylor (2002); e: Weisberg & Taylor (2004); f: Stairs et al. (2002); 
g: Bogdanovet al. (2002); h: Stairs. Thorsett & Arzoumanian (2004) 
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Previously, the moment of inertia of the Crab pulsar was estimated by Bejger & Haensel (2003). They 
used an estimate (Fesen, Shull & Hurford 1997) for the mass of the ionized portion of the Crab's remnant, 
4.6 ±1.8 M Q , to infer a lower limit to the Crab pulsar's moment of inertia of 97 ±38 M Q km 2 . Within errors, 
the lower limit rules out only the softest equations of state, as is shown in Fig. 1. 

It should be noted that the binary PSR J0737-3039, although highly relativistic, is not favorably inclined 
for observation of the precession of the inclination, which is proportional to cosi. But this particular 
inclination does allow a simplification of the interpretation of the pcriastron advance. Not only is the spin- 
orbit contribution nearly independent of the pulsar orientation angle 4>a, but errors in i are also not going 
to be a significant restriction to a measurement of I a- 

It is interesting to compare PSR 0737-3039 with other known relativistic binaries. Table 1 compares 
properties of the relativistic binaries for which spin orientation of pulsar A can be estimated. It is noteworthy 
that the net timing delays caused by precession-induced inclination shifts are more than an order of magnitude 
larger for PSR 1913+16 and PSR 1534+12 than for PSR 0737-3039, due to their smaller inclinations. 
However, the precessional periods of these two systems are 4 times larger, and their periastron advances are 
about 6 times smaller, than for PSR 0737-3039. This factor of 24 in observability is significant, and explains 
why measurements of this effect in these systems have not been made. 

The nearness and faintness of this binary gives hope that other highly relativistic systems might even- 
tually be observed. If further highly relativistic systems are discovered, it seems unlikely that any of them 
will have an inclination angle so unfortunate for detecting the precession of the orbital plane. With two 
spin-orbit effects to be observed, the moment of inertia might be measured to even greater accuracy than 
contemplated for PSR 0737-3039. 
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